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ABSTRACT: The recent literature shows a renewed interest, with various independent 
approaches, in the classical models for spin. Considering the possible interest of those 
results, at least for the electron case, we purpose in this paper to explore their physical 
and mathematical meaning, by the natural and powerful language of Clifford algebras 
(which, incidentally, will allow us to unify those different approaches). In such models, the 
ordinary electron is in general associated to the mean motion of a point-like "constituent" 
Q, whose trajectory is a cylindrical helix. We find, in particular, that the object Q obeys 
a new, non-linear Dirac-like equation, such that — when averaging over an internal cycle 
(which corresponds to a linearization) — it transforms into the ordinary Dirac equation 
(valid, of course, for the electron whole) . 
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1. Introduction - The possibility of constructing formal classical models of spin was 
already realized at least 60 years ago, from different points of view. 1 In particular, 
Schrodinger's suggestion 2 that the electron spin was related to its Zitterbewegung (zbw) 
motion has been investigated by several authors 3 . 

In ref. 4 , for instance, one meets even the proposal of models with clockwise and anti- 
clockwise "inner motions" as classical analogues of quantum relativistic spinning particles 
and antiparticles, respectively. The use of grassmannian variables in a classical lagrangian 
formulation for spinning particles was proposed, on the contrary, by Berezin and Marinov 5 
and by Casalbuoni 6 . Pauri, 7 moreover, showed how a relativistic (or non-relativistic) spin 
can enter a classical formulation, in the ordinary phase space [without resource now to 
grassmannian quantities], just as a consequence of the algebraic structure of the Poincare 
(or Galilei) group. In such an interesting approach, it was found 7 that a relativistic classi- 
cal zbw motion can directly follow as a "spin effect" from the requirements of existence of 
a lagrangian and of a covariant position in Dirac's instant form of relativistic dynamics. 
The quantum analogue of those developments had been studied in ref. 8 . 

The number of papers appeared on the subject of classical models for spin, starting 
from the fifties, is so large that it would be difficult to try quoting them here. A recent 
approach, based on a generalization of Dirac non-linear electrodynamics, where the spin of 
the electron is identified with the momentum of the Poynting vector in a soliton solution 
of that theory, can be found in ref. 9 . 

In this paper we choose, by making recourse to the natural and powerful language of 
Clifford algebras, to refer ourselves mainly to the model by Barut and Zanghi 10 ' 11 (BZ), 
which relates the spin (at least in the case of the electron) to a helical motion. Namely, 
we first recast the BZ model into the Clifford formalism, in the meantime clarifying its 
physical and mathematical meanings; then, we quantize that model for the electron case. 
In particular, we derive from the BZ model a new non-linear equation for the "spinorial 
variable" of the model, which when linearized reduces to the ordinary Dirac equation. 
Solutions of this non-linear equation will be discussed in forthcoming papers 12 . 



2. Spin and helical motion - In the BZ model, 10 a classical electron is characterized, 
besides by the usual pair of conjugate variables (x M ,p M ), by a second pair of conjugate 
classical spinor variables (z,iz), representing internal degrees of freedom, which are func- 
tions of an invariant time parameter r, that when convenient will be identified with the 
proper time of the center of mass (CM). Quantity z was a Dirac spinor, while z = z'~f°. 
Barut and Zanghi, then, introduced for a spinning particle the classical lagrangian [c = 1] 

C = \\i{zz — zz) + p^x^ — z^z) + eA fl (x)zj^z , (1) 

where A has the dimension of an action, 7 M are the Dirac matrices, and the particle 
velocity is 
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Vp^zi^z . (V) 

We are not writing down explicitly the spinorial indices of z and z. Let us consider 
the simple case of a free electron (A^ = 0). Then a possible solution of the equations of 
motion corresponding to the lagrangian (1) is: #1 



v 

z(t) = [cos mr — 27^ — sin mr] z(0) , (2a) 

m 

V 

z(t) = z(0) [cos mr + 27^— sin mr] , (2b) 



m 



and p M = constant; p = m ; H = p^zj^z = p^v M ; and finally: 



± = v = EJLH + \xJ0) - -^#]cos2mr + ^-Uin2mr (2c) 
m 2 m 2 2m 



(which in ref. 10 appeared with two misprints). In connection with this "free" general 
solution, let us remark that if is a constant of motion so that we can set H = m. 
Solution (2) exhibits the classical analog of the phenomenon known as "zitterbewegung" 
(zbw): in fact, the velocity t> M = x^ contains the (expected) term p^/m plus a term 
describing an oscillatory motion with the characteristic frequency uo = 2m. The velocity 
of the center of mass will be given by = p^/m. Notice incidentally that, instead of 
adopting the variables z and z, one can work in terms of the spin variables, i.e. in terms 
of the dynamical variables (x^, v^, 71^, S^), where: 
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vr M = Pf, - eA^ = -izYi^ ~{ u }z , 



so that S^v = ft pV v — tt u v^ and = AS^ir". In the case of a free electron, by varying 
the action corresponding to C one finds as generator of space-time rotations the conserved 
quantity J^ v = L^ v + S^ u , where S^ u is just the particle spin.* 2 

Let us explicitly observe that solution (2c) is the equation of a space-time cylindrical 
helix, i.e. it represents in 3-space a helical motion. Let us also stress that this motion 
describes particle spin at a classical level. In fact, such a classical system has been shown 
by Barut and Pavsic 13 to describe, after quantization, the Dirac electron. 13 ~ 15 

An alternative approach leading to a classical description of particles with spin was 
put forth by Pavsic, 16-18 by making recourse to the (extrinsic) curvature of the particle 
world-line in Minkowski space. Let us finally mention that the same classical equations of 

o#i p or other solutions, see refs. 12 . 

°# 2 Alternatively, Barut and Zanghi, 10 in order to study the internal dynamics of the considered 
(classical) particle, did split x^ and = as follows: x^ = + Q^; w M = + (where by 
definition = and = Q^). In the particular case of a free particle, — 0; W M — p^/m. One 
can now interpret X^ and p M as the cm. coordinates, and and P M = ml) ' ^ as the relative position and 
momentum, respectively. For a free particle, then, one finds that the internal variables are coordinates 
oscillating with the zbw frequency 2m; and that, again, the total angular momentum J^ v = L M „ + is 
a constant of motion, quantities S^ v being the spin variables. 
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motion (and the same Poisson-bracket algebra) have been found also in a third approach, 
which consists in adding to the ordinary lagrangian an extra term containing Grassmann 
variables. 19 



3. About the electron structure - Considering the interest of the previous results (which 
suggest in particular that the helical motion can have a role in the origin of spin), we pur- 
pose to explore their physical meaning more deeply, by the very natural — and powerful- 
language of the Clifford algebras: 20 ' 21 in particular of the "space-time algebra (STA)" 
Mifi. First of all, let us preliminarily clarify why Barut and Zanghi had to introduce the 
Dirac spinors z in their lagrangian, by recalling that classically the motion of a spinning 
top has to be individuated by (i) the world-line a of its center of mass {e.g., by the 
ordinary coordinates x^ and the conjugate momenta and (ii) a Frenet tetrad* 3 
attached 22 to the world-line a. This continues to be true when wishing to describe the 
motion of a point-like spinning particle. For the Frenet tetrad 23 we have: 

e, = R llx R = A£ 7v ; A;gL| (3) 

where e is parallel to the particle velocity v (even more, e = v whenever one does use 
as parameter r the CM system proper-time); the tilde represents the reversion^ 4 ; and 
R = R(t) is a "Lorentz rotation" [more precisely, R G Spin + (1,3), and a Lorentz 
transform of quantity a is given by a' = RaR] . Moreover RR = RR = 1 . The Clifford 
STA fundamental unit-vectors 7 M should not be confused with the Dirac matrices 7 M . Let 
us also recall that, while the orthonormal vectors 7 M = d/dx^ constitute a global tetrad 
in Minkowski space-time (associated with a given inertial observer), on the contrary the 
Frenet tetrad e M is defined only along a, in such a way that e is tangent to a. At last, 
it is: 7^ = rf v ^ v , and 75 = 7o7i7273- 

Notice that R(t) does contain all the essential information carried by a Dirac spinor. 
In fact, out of R, a "Dirac-Hestenes" (DH) spinor 24 -^dh can be constructed as follows: 

^dh = P h^/ 2 R (4) 

where p is a normalization factor; and e' 375 = +1 for the electron (and —1 for the 
positron); while, if e is a primitive idempotent of the STA, any Dirac spinor ipD can be 
represented in our STA as: 25 

V'd = ^dh£. (5) 



o#3 r-Q^ uge £ p rene ^ tetrads in connection with the Dirac formalism was first investigated in ref. 22 . 
o#4 r-Q^ mam an ti-automorphism in M1.3 (called reversion), denoted by the tilde, is such that 
AB = BA, and A = A when A is a scalar or a vector, while F — —F when F is a 2-vector. 
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For instance, the Dirac spinor z introduced by BZ is obtained from the DH spinor* 5 by the 
choice e = 1(1 + 70). Incidentally, the Frenet frame can also write pe M = ^dhT^dh- 

Let us stress that, to specify how does the Frenet tetrad rotate as r varies, one has 
to single out a particular R(r), and therefore a DH spinor -?/> D h , and eventually a Dirac 
spinor ip D . This makes intuitively clear why the BZ Dirac-spinor z provides a good 
description of the "spin motion" of a classical particle. 

Let us now repeat what precedes on a more formal ground. In the following, unless 
differently stated, we shall indicate the DH spinors -^dh simply by ip. 



Let us translate the BZ lagrangian into the Clifford language. 

In eq.(l) quantity z T = (z\ Z2 z 3 z±) is a Dirac spinor r — > z(t), and z = 2^7°. 
To perform our translation, we need a matrix representation of the Clifford space-time 
algebra; this can be implemented by representing the fundamental Clifford vectors (70, 
7i? 72, 73) by the ordinary Dirac matrices 7 . Choosing: 




7o >7o = i n n I : 7 — - 7, 

the representative in Mi^ of Barut-Zanghi's quantity z is 

z ^^ = ^ e] e =i(i + 7o ) 

where ip, and ip, are represented (with ipip = 1) by 



1> 



The translation of the various terms in eq.(l) is then: 




( Z1 


-h 


Z 3 


Z~4 


\ 




( 


Zl 


Z~2 


-Z3 


-Z4 \ 


Z2 


Zl 


Z4 


-Z 3 








-z 2 


Zl 


-Z4 


Z3 




Z~4 


Zl 


-z 2 






-Z3 


-Z4 


Zl 


-z 2 


V Zi 


-Z 3 


Z2 


Zl 


J 




{ 


-Z4 


Z3 


-Zl 


Zl ) 



\i{zz 



(6) 



(7) 



zz) — > (^^7i72)o 
- zYz) — >• (p(x - 1pjQlp)) 
eA^zYz — ► e (A^7o^)o , 

where ( ) means "the scalar part" of the Clifford product. Thus, the lagrangian £ in 
the Clifford formalism is 



C = (?/n/>7i72 + p(x - ^70^) + eA^7o^)o , (8) 

which is analogous, incidentally, to Kriiger's lagrangian 26 (apart from a misprint). 

As we are going to see, by "quantizing" it, also in the present formalism it is possible 

o#5 rpjjg pjH spinors can be regarded as the parent spinors, since all the other spinors of common use 
among physicists are got from them by operating as in eq.(5). We might call them "the fundamental 
spinors" . 
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(and, actually, quite easy) to derive from C the Dirac-Hestenes equation:* 6 



dtj}{x) 7172 + mip(x)^Q + eA(x)ip(x) = 0, (9) 

which is nothing but the ordinary Dirac equation written down in the Clifford 
formalism. 20 ' 24 Quantity d = 7 M <9 M is the Dirac operator. Let us notice that p in 
eq.(8) can be regarded as a Lagrange multiplier, when the velocity v — x is represented 
by ip'jo'ip. The BZ model is indeed a hamiltonian system, as proved by using Clifford 
algebras in ref. 27 (cf. also ref. 28 ). The dynamical variables are then (■0, ip, x,p), and the 
Euler-Lagrange equations yield a system of three independent equations: 

■07i 72 + 7T07o _= (10a) 
x = ■07oV' (10b) 
■k = eF-x (10c) 

where F = d A A is the electromagnetic field (a bivector, in Hestenes' language) and 
7i = p — eA is the kinetic momentum. [Notice incidentally, from eq.(lOb), that x 2 = p(r)]. 

At this point, let us consider a velocity vector field V(x) together with its integral lines 
(or stream-lines). Be a the stream-line along which a particle moves (i.e., the particle 
world-line). Then, the velocity distribution V is required to be such that its restriction 
V(x): to the world-line a is the ordinary velocity v = v(r) of the considered particle. 

If we moreover recall 29 ' 30 that any Lorentz "rotation" R can be written R = e^, where 
T is a bivector, then along any stream-line a we shall have: 19 

R=^l = - v »n„R = -SIR, (11) 
dr 2 M 2 v ; 

with d^R = VL^R/2, where f2 M = Id^T , and where il = is the angular-velocity 

bivector (also known, in differential geometry, as the "Darboux bivector"). Therefore, 
for the tangent vector along any line a we obtain the relevant relation: 

= v^dn — v ■ d 

dr 

The [total derivative] equation (10a) thus becomes:* 6 

v ■ (9-07i72 + vr-070 = , (12) 

which is a non-linear [partial derivative] equation, as it is easily seen by using eq.(10b) 
and rewriting it in the noticeable form 



°# 6 Observe that in eq.(8) it is ip = tp(r), while in eq.(9) we have ip — ip(x) with ip(x) such that its 
restriction VK x )|cr t° the world-line a coincides with iP(t). Below, we shall meet the same situation, 
for instance, when passing from eq.(10a) to eqs.(12)-(12'). 
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(V'ToV') • <9^7i72 + vr^7o = 



(12') 



Equation (12') constitutes a new non-linear Dirac-like equation. The solutions of this 
new equation will be explicitly discussed in refs. 12 . Let us here observe only that the 
probability current J = V is conserved: d • V = 0, as we shall show elsewhere. 
Let us pass now to the free case (A^ = 0), when eq.(lOa) may be written 

^7i72 + P^lo = , (10'a) 

and admits some simple solutions. 12 Actually, in this case p is constant [cf. eq.(lOc)] and 
one can choose the 7 M frame so that p = m7 is a constant vector in the direction 7 . 
Since x = V'ToV'j it follows that 

1 ~ p 

v = — ■ippip ; — = ifj^vip^ 1 . (13) 
m m 

The mean value of v over a zitterbewegung period is then given by the relation 

(v) z bw = p/m = ^)- l vi)~ l (14) 

which resembles the ordinary quantum-mechanical mean value for the wave-function -?/> _1 
(recall that in Clifford algebra for any ip it exists its inverse). Let us recall, by comparison, 
that the time average of given in eq. (2c) over a zbw period is evidently equal to p^j m. 

Let us explicitly stress that, due to the first one [z — > ipe] of eqs.(6), the results 
found by BZ for z are valid as well for if) in our formalism. For instance, for BZ [cf. 
eq.(l')] it was = zj^z = v£ z , while in the Clifford formalism [cf. eq.(7)] it is 
v = = (7°^7 M ^)o7^ = v Bzlfi- ^ s a consequence, a refers in general to a 

cylindrical helix (for the free case) also in our formalism. 

Going back to eq.(lO'a), by the second one of eqs.(13) we finally obtain our non-linear 
(free) Dirac-like equation in the following form: 

v ■ <9?/>7i72 + mip^vip^ip^Q = , (15) 

which in the ordinary, tensorial language would write: i(^/^f tJ '^/)d^ = "y^p^, where 
^ and 7 M are now an ordinary Dirac spinor and the ordinary Dirac matrices, respectively, 
and p^ = id ^ 

In connection with this fundamental equation of motion (15), let us explicitly notice the 
following. At a classical level, the equation of motion in the BZ model was eq.(lOa), which 
held for the world-line a. In other words, eq.(lOa) was valid for one world-line; on the 
contrary, eq.(15) is a field equation, satisfied by quantities ip(x) such that ip(x)^ a = ip(r). 
A change in interpretation is of course necessary when passing from the classical to the 
"quantum" level: and therefore eq.(15) is now to be regarded as valid for a congruence of 
world lines, that is to say, for a congruence of stream-lines of the velocity field V = V(x). 
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In the quantum case, the "particle" can follow any of those integral lines, with probability 
amplitude p. In this context, it must be recalled that a tentative interpretation of the 
Dirac equation within the Clifford algebra approach has been suggested in ref. 20 , and 
later in ref. 31 . However, in the present paper we shall not put forth, nor discuss, any 
interpretation of our formalism. 

As we have seen, eq.(15) will hold for our helical motions. Notice moreover that, since 
[cf. eq.(4)] it is ip — p^e^^R, in the case of the helical path solution the Lorentz 
"rotation" R will be the product of a pure space rotation and a boost. 

It is important to observe that, if we replace v in eq.(15) by its mean value over a zbw 
period, eq.(14), then we end up with the Dirac-Hestenes equation [i.e., the ordinary Dirac 
equation!], valid now for the center-of-mass world-line. In fact, since (f) z b w = p/m, 
eq.(15) yields 

p ■ <9^7i72 + = , (15') 

and therefore — if we recall that for the eigenfunctions of p in the DH approach 20 it holds 
•9^7172 = VPi so that (p ■ (9)^7i72 = P<9^7i72 — one obtains: 

p(<9^7i7 2 + m^7o) = 0, (15") 

which is satisfied once it holds the ordinary Dirac equation (in its Dirac-Hestenes form): 

•9^7172 + «#7o = . (16) 

Let us observe that all the eigenfunctions of p are solutions both of eq.(15) and of the 
Dirac equation. 

In conclusion, our non-linear Dirac-like equation (15) is a quantum-relativistic equa- 
tion, that can be regarded as "sub-microscopic" in the sense that it refers to the internal 
motion of a point-like "constituent" Q. In fact, the density current ip^f°ip, relative to the 
solutions ifj of eq.(15), does oscillates in time in a helical fashion, in complete analogy 
to the initial equation (2c); so that x and x refer to Q. On the contrary, the ordinary 
(linear) Dirac equation can be regarded as the equation describing the global motion 
of the geometrical centre of the system (i.e., of the whole "electron"); actually, it has 
been obtained from eq.(15) by linearization, that is to say, replacing the density current 
v = ip'-f ^ by its time average p/m over the zbw period [cf. eqs.(13)-(14)]. 

At last, let us underline that, in the free case, eq.(lO'a) admits also a trivial solution 
(j , corresponding to rectilinear motion. 



4. A very simple solution - In the free case, the equation (10'a) admits also a very 
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simple solution (the limit of the ordinary helical paths when their radius r tends to zero), 
which — incidentally — escaped BZ's attention; namely: 

tp = p2 exp [—7271 mr] = ijj(0) exp[— , j2li mr ] (17) 

as it can be easily verified. Notice that expressions (17) are also solutions of the Dirac 
equation (16). Moreover, from eq.(lOb) it follows that v = x = p7o and we can set 
p = 1, which confirms that the trivial solution (17) corresponds to rectilinear uniform 
motion, and that r in this case is just the proper-time along the particle world-line &q. 
(In this case, of course, the zbw disappears). 

Before going on, we want to explicitly put forth the following observation. Let us first 
recall that in our formalism the (Lorenz force) equation of motion for a charged particle 
moving with velocity w in an electromagnetic field F is 

w = —F-w. (18) 
m 

Now, for all the free-particle solution of the BZ model in the Clifford language, it holds 
the "Darboux relation" : 

e M = Q-e^, (19) 

so that the "sub-microscopic" point-like object Q, moving along the helical path a, is 
endowed [cf. eq.(ll)] with the angular-velocity bivector 

n = \ e M A = \ e„e" , (20) 

as it follows by recalling that e M can always be written, like in eq.(3), as e M = R^f^R. 
Finally, let us observe that eq.(19) yields in particular eo = Q ■ eo, which is formally 
identical to eq.(18). Thus, the formal, algebraic way we chose [see eq.(18)] for describing 
that the system as a whole possesses a non-vanishing magnetic dipole structure suggests 
that the bivector field Q may be regarded as a kind of internal electromagnetic-like field, 
which keeps the "sub-microscopic" object Q moving along the helix. 32-36 In other words, 
Q may be considered as confining itself along a [i.e., along a circular orbit, in the electron 
CM.], via the generation of the internal, electromagnetic-like field 

Tfl 

F int = — e„Ae M . 
We shall further discuss this point elsewhere. 



5. Further remarks - In connection with our new equation (12'), or rather with its 
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(free) form (15), we met solutions corresponding — in the free case — to helical motions 
with constant radius r; as well as a limiting solution, eq.(17), for r — > 0. We have seen 
above that the latter is a solution also of the ordinary (free) Dirac equation. Actually, 
the solution of the Dirac equation for a free electron in its rest frame can be written 20 in 
the present formalism as if)(x) = i/)(0) exp[— 7271 mr], which coincides with our eq.( 17) 
along the world-line &q. Let us here add the observation that, by recalling that in the 

Clifford formalism the spin bivector S is given by S = -R7271-R— = e 2 ei — , whilst the 
angular-velocity bivector Q is given by eq.(20), one can finally write 

p. v = n-S = m, (21) 

which seems to suggest 20 the electron rest-mass to have an (internal) kinetic origin! It is 
possible, moreover, that such a motion be also at the origin of electric charge itself; in any 
case, we already saw that, if the electron is associated with a clock-wise rotation, then the 
positron will be associated with an anti-clock-wise rotation, with respect to the motion 
direction. Further considerations about the solutions of our non-linear, new equation and 
the interesting consequences of the present formalism are in preparation and will appear 
elsewhere. 

At last, we mentioned at the beginning of this note about further methods for intro- 
ducing an helical motion as the classical limit of the "spin motion" . We want here to 
show how to represent in the Clifford formalism the extrinsic curvature approach, 16 ' 17 due 
to its possible interest for the development of the present work. 

Let us first recall that in classical differential geometry one defines 23 the Frenet frame 
{e M } of a non-null curve a by the so-called Frenet equations, which with respect to proper 
time r write [besides x — eo — v\: 

x = e = K 1 e 1 ; h = -K.e + K 2 e 2 ; e 2 = -K 2 e l + K 3 e 3 ; e 3 = -K 3 e 2 (22) 

where the i-th curvatures K x (i=l,2,3) are scalar functions chosen in such a way that 
ej 2 = —1, with j = l,2,3. Quantity K\ is often called curvature, and K 2 , K3 torsions 
(recall that in the 3-dimensional space one meets only K\ and K 2 , called curvature and 
torsion, respectively). Inserting eqs.(22) into eq.(20), we get for the Darboux (angular- 
velocity) bivector: 

Q = XieV + K 2 e 2 e l + K 3 e 3 e 2 , (23) 

so that one can build the following scalar function 

Q-Q = K\-K\-Kl = (e M A e M ) • (e u A e u ) . (23') 

At this point, one may notice that the square, K 2 , of the "extrinsic curvature" entering 
Pavsic's lagrangian is equal to —K 2 , so that such a lagrangian 16 results — after the present 
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analysis — to take advantage only of the first part, 

x 2 = el = -K\ , 

of the Lorentz invariant (23'). On the contrary, in our formalism the whole invariant fl-fl 
suggests itself as the suitable, complete lagrangian for the problem at issue; and in future 
work we shall exploit it, in particular comparing the expected results with Plyushchay's. 37 
For the moment, let us stress here only the possibly important result that the lagrangian 
£ = Q ■ Q does coincide (factors apart) along the particle world-line a with the auto- 
interaction term 38 

e 5 (d^ a M ) • (dr a 0„) 

of the Einstein-Hilbert lagrangian density written (in the Clifford bundle formalism) in 
terms of tetrads of 1-form fields 9^. Quantity 9 5 = 9°9 1 9 2 9 3 is the volume element. 

Finally, we can examine within our formalism the third approach: that one utilizing 
Grassmann variables. 19 For instance, if we recall that the Grassmann product is nothing 
but the external part A T A B s = (A T B s )\ r _ s \ of the Clifford product (where A T , B s are a 
r- vector and a s- vector, respectively), then the Ikemori lagrangian 19 can be immediately 
translated into the Clifford language and shown to be equivalent to the BZ lagrangian, 
apart from the constraint p 2 = m 2 . 
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